Mon. Not. R. Astron. Soc. 000, 000-000 (0000) 



Printed 29 August 2008 



(MN M£X style file v2.2) 



Exploring the substellar temperature regime down to 
-550K 



oo ! Ben Burningham 1 *, D. J. Pinfield 1 , S. K. Leggett 2 , M. Tamura 3 , P. W. Lucas 1 , 
§ ;D. Homeier 4 , A. Day- Jones 1 , H. R. A. Jones 1 , J.R.A. Clarke 1 , 

; M. Ishii 5 , M. Kuzuhara 6 , N. Lodieu 7 , M. R. Zapatero Osorio 7 , B.P. Venemans 8 , 
2f: D.J. Mortlock 9 , D. Barrado y Navascues 10 ,E. L. Martin 7 , A. Magazzu 11 

1 Centre for Astrophysics Research, Science and Technology Research Institute, University of Hertfordshire, Hatfield ALIO 9AB 

2 Gemini Observatory, 670 N. A 'ohoku Place, Hilo, HI 96720, USA 
3 National Astronomical Observatory, Mitaka, Tokyo 181-8588 

4 Institut fur Astrophysik, Georg- August- Universitat, Friedrich-Hund-Platz 1, 37077 Gottingen, Germany 

5 Subaru Telescope, 650 North A 'ohoku Place, Hilo, Hi 96720, USA 

6 University of Tokyo, Hongo, Tokyo 113-0033, Japan 

i 7 Instituto de Astrofisica de Canarias, 38200 La Laguna, Spain 
O i 8 Institute of Astronomy, Madingley Road, Cambridge CB3 0HA, UK 

9 Astrophysics Group, Imperial College London, Blackett Laboratory, Prince Consort Road, London SW7 2AZ 
O 1 10 Laboratorio de Astrofisica Espacial y Fisica Fundamental, INTA, P.O. Box 78, E-28691 Villanueva de la Canada (Madrid), Spain 
11 Fundacion Galileo Galilei-INAF, Apartado 565, E-38700 Santa Cruz de La Palma, Spain 



3 

(N 



C3 



(N 
> 

\o 
o 
o 

o 

oo 
o 



1—1 



29 August 2008 



ABSTRACT 

We report the discovery of three very late T dwarfs in the UKIRT Infrared Deep 
Sky Survey (UKIDSS) Third Data Release: ULAS J101721.40+011817.9 (ULAS1017), 
ULAS J123828.51+095351.3 (ULAS1238) and ULAS J133553.45+113005.2 
(ULAS1335). We detail optical and near-infrared photometry for all three sources, and 
mid-infrared photometry for ULAS1335. We use near-infrared spectra of each source 
to assign spectral types T8p (ULAS1017), T8.5 (ULAS1238) and T9 (ULAS1335) 
to these objects. ULAS1017 is classed as a peculiar T8 (T8p) due to appearing as 
a T8 dwarf in the J-band, whilst exhibiting H and A-band flux ratios consistent 
with a T6 classification. Through comparison to BT-Scttl model spectra we estimate 
that ULAS1017 has 750K<T eff <850K, and 5.0 < log g(cm S - 2 ) < 5.5, assuming 
solar metallicity. This estimate for gravity is degenerate with varying metallicity. 
We estimate that ULAS1017 has an age of 1.6-15 Gyr, a mass of 33-70 Mj and 
lies at a distance of 31-54 pc. We do not estimate atmospheric parameters for 
ULAS1238 due to a l ack of A-band photometry. We extend the unified scheme of 
iBurgasser et al.l ([2006) to the type T9 and suggest the inclusion of the Wj index to 
replace the now saturated J-band indices. We propose ULAS1335 as the T9 spectral 
type standard. ULAS1335 is the same spectral type as ULAS J003402. 77-005206.7 
and CFBDS J005910. 90-011401. 3. We argue that given the similarity of the currently 
known >T8 dwarfs to the rest of the T dwarf sequence, the suggestion of the Y0 
spectral class for these objects is premature. Comparison of model spectra with 
that of ULAS1335 suggest a temperature below 600K, possibly combined with 
low-gravity and/or high- metallicity. We find ULAS 1335 to be extremely red in 
near to mid-infrared colours, with H — [4.49] = 4.34 ± 0.04 . This is the reddest 
near to mid-infrared colour yet observed for a T dwarf. The near to mid-infrared 
spectral energy distribution of ULAS1335 further supports T e g < 600K, and we 
estimate T cS ~ 550 - 600K for ULAS1335. We estimate that ULAS1335 has an age 
of 0.6-5.3 Gyr, a mass of 15-31 Mj and lies at a distance of 8-12 pc. 
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1 INTRODUCTION 



2.1 Mining UKIDSS DR3 



Extending the known sample of field dwarfs to ever lower 
effective temperatures (T e ff) is important not only for the 
determination of the field mass function, but also for prob- 
ing temperature and pressure regimes that have hitherto 
been unexplored observationally. The study of extremely 
cool brown dwarfs opens a window on atmospheric physics 
which will be fundamental for understanding the processes 
within a broad range of substellar atmospheres, including 
those of giant exoplanets. 

At the time of writing there are 142 T d warfs published 
in th e literature (e.g., DwarfArchives.org; IPinfield et al.l 
2008). Of these, onl y five are class i fied as T8 or later 
(using the scheme of iBurgasser et all [2006). two of which 
have spectral types later than T8: ULAS J003 4 02.77- 
005206.7 (hereafter ULAS0034; IWarren et all 120071 ) 
and CFBDS J0059 10.90-011401.3 (hereafter CFBDS0059; 
iDelorme et"al]|2008l). ULA S0034 has an inferred T eff o f 600- 
700 K |Warren et al.ll2007l ). and lDelorme et all (|2008D infer 
a T off for CFBDS0059 that is ~ 50 K cooler than that of 
ULAS0034. 

Expanding the sample of very cool brown dwarfs, 
and exploring possible factors that could motivate the 
implementation of new spectral classes (e.g Y dwarfs; 
iKirkpatrick et alj Il999l) are primary science drivers for 
the UKIRT Infrared De ep Sky Survey (UKID SS) Large 
Area Survey (LAS; see lLawrence et al] 120071 ). Our col- 
laboration, the UKIDSS Cool Dwarf Science Working 
Group (CDSWG) is e ngaged in a substantial effort to 
achieve this aim (e. g.. IPinfield et all Ejoll; iKendall et al.l 
l2007l ; lLodieu et all I2007I ; I Warren et al] 120071 ). Here we 



present the recent discovery by CDSWG of three very 
cool brown dwarfs. These are ULAS J101721. 40+011817.9 
(ULAS1017), ULAS J123828.51+095351.3 (ULAS1238) and 
ULAS J133553.45+113005.2 (ULAS1335). The last of these, 
ULAS1335, may be the coolest brown dwarf yet discovered 
and we here discuss the properties of this object, estimating 
a T e ff of approximately 550-600 K. 



2 SOURCE IDENTIFICATION 

The objects presented here were selected from the third data 
release (DR3) of UKIDSS (Warren et al., in prep). We used 
the same method as that outlined in Sections 2 and 3 of 
IPinfield et all (2008). The process of identifying and con- 
firming extremely late T dwarfs can be separated into three 
principle stages: 1) initial selection by mining survey data; 
2) confirmation of photometric properties by near-infrared 
and optical follow-up; 3) near-infrared spectroscopy to esti- 
mate spectral type. The Gemini/NIRI spectroscopic follow- 
up took place in two stages in queue mode. Initial, short 
(tint = 16mins) J-band observations were obtained first. 
These initial spectra confirmed that two of our targets were 
indeed of extremely late spectral type, and deeper J, H and 
if-band observations were then queued. Subaru/IRCS spec- 
troscopic follow-up took place separately for additional can- 
didates in classical mode and confirmed the very late type of 
a third target. We will discuss each of these stages in more 
detail in the following sections. 



We initially searched the UKIDSS LAS data for sources that 
met colour criteria: Y — J > 0.5; J — H < 0.1; a < 0.3, where 
a is the uncertainty in the magnitude measured in each 
of the three yj_ff-bands. The resultant candidate list was 
cross-matched against the Sloan Digital Sky Survey (SDSS) 
Data Release 6 (DR6) using a pairing radius of 2 arcsecs. 
Objects with no SDSS counterpart, or whose SDSS-UKIDSS 
z — J > 3.0, were retained as potential candidates. We iden- 
tified some 50 candidate T dwarfs in this search. ULASJ1238 
and ULASJ1335 were picked out for immediate and rapid 
follow-up due to survey colours that were reminiscent of the 
prototype >T8 dwarf ULAS J0034. ULAS1017 on the other 
hand was selected due to its extremely blue J — H colour 
from UKIDSS photometry. 



2.2 Ground based photometric follow-up 

Near infrared follow-up photometry was obtained usin 
the UKIRT Fast Track Imager (UFTI; iRoche et all 1200; 
mounted on UKIRT, and the Long-slit Infr ared Imaging 
Spectrograph (LIRIS; iManchado et al.l ll998T ) mounted on 
the William Herschel Telescope on La Palma. The mosaics 
were produced using sets of jittered images, with individual 
exposure times, jitter patterns and number of repeats given 
in Table [T] The data were dark subtracted, Hatfield cor- 
rected, sky subtracted and mosaiced using ORAC-DR for 
the UFTI data, and LIRIS-DR for the LIRIS data. We cali- 
brated our UFTI obs ervations using UKIRT Faint Standards 
l|Leggett et al] |2006). with a standard observed at a similar 
airmass for each target. All UFTI data were obtained under 
photometric conditions, with seeing better than 0.9 arcsec. 
The wider field LIRIS data were obtained through thin cir- 
rus in variable seeing (0.8 - 2.0 arcsec) and were calibrated 
against 2MASS stars within the field of view. LIRIS uses 
a K s filter, and as such no transform was required for the 
2MASS stars. The K s magnitude for ULAS1017 was trans- 
formed to the MKO system using the transform derived by 
IPinfield et ail (|2008T ). 

We have also obtained optical z-band photometry using 
the ESO Multi-Mode Instrument (EMMI) mounted on the 
New Technology Telescope at La Silla, Chile on the nights of 
2008 January 29 and 2008 January 30. These observations 
are summarised in Table[T] For this optical follow-up we used 
a Bessel z-band filter (ESO Z#611). The data were reduced 
using standard IRAF packages, and then multiple images of 
the same target were aligned and stacked to increase signal- 
to-noise. 

We calibrated our z— band photometry using SDSS DR6 
sources present in our images as secondary standards. We 
transformed the Sloan z'(AB) magnitu des to the EMMI sys - 
tem using the transformation given by I Warren" et all (|2007l) 
before using them to determine zero points, which had a typ- 
ical scatter of ±0.04 mags. The resulting EMMI photometry 
for our tar gets was then trans formed to the Sloan z'(AB) 
system (see Warren et al]|2007l ). The results of our ground- 
based follow-up photometry are given in Table [2] which also 
gives the original database photometry. 
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Object 


Filter 


Instrument 


UT 


Date 


Total integration time 


tint 


breakdown 


ULAS1017 


,/ 


UFTI 


2008 


Jan 17 


300s 


Ci — 5 r 


— 1 ; ^exp 


= 60s) 




H 


UFTI 


2008 


Jan 17 


2400s 


(j = 5, r 


— 6, iexp 


= 60s) 




K s 


LIRIS 


2008 


Mar 17 


2400s 


(j = 5, r = 


- 24, £ cx p 


= 20s) 




Z EMMI 


EMMI 


2008 


Jan 30 


2400s 


(j = l,r = 


4, £exp - 


= 600s) 


ULAS1238 


J 


UFTI 


2008 


Jan 25 


300s 


(j = 5, r 


— 1; ^CXp 


= 60s) 




H 


UFTI 


2008 


Jan 25 


1800s 


(j = 5, r 


— 6, iexp 


= 60s) 




ZEMMI 


EMMI 


2008 


Jan 30 


2400s 


(j = V = 


4, iexp - 


= 600s) 


ULAS1335 


Y 


UFTI 


2008 


Jan 16 


540s 


(j = 9, r 


— 1 ; iexp 


= 60s) 




J 


UFTI 


2008 


Jan 16 


300s 


(j = 5, r 


— 1 , texp 


= 60s) 




H 


UFTI 


2008 


Jan 16 


900s 


(j = 5, r 


— 3, f-exp 


= 60s) 




K 


UFTI 


2008 


Jan 16 


900s 


(j = 5, r 


— 3, f. e xp 


= 60s) 




Z EMMI 


EMMI 


2008 


Jan 29 


1200s 


(j = l,r = 


2, £exp - 


= 600s) 



Table 1. Summary of the observations obtained for near infrared and optical photometric follow-up. The breakdown of each integration 
is given in the final column with the following notation: j = number of jitter points; r = number of repeats for jitter pattern; t CX p = 
exposure time at each jitter point. 



Object Y J H K z'(AB) z'(AB)-J Y-J J-H H-K 



ULAS1017 


19.72 


± 


0.15 


18.57 


± 


0.07 


19.38 


± 


0.24 




















1.15 


± 


0.17 


-0.82 ±0.25 












18.53 


± 


0.02 


19.07 


± 


0.02 


19.01 


± 


0.10 


22.49 


± 


0.09 


3.96 


± 


0.09 


1.19 


± 


0.15 


-0.54 ±0.03 


0.06 ±0.10 


ULAS1238 


19.55 


± 


0.07 


18.78 


± 


0.06 


19.25 


± 


0.15 




















0.77 


± 


0.09 


-0.47 ±0.17 












18.95 


± 


0.02 


19.20 


± 


0.02 








22.74 


± 


0.12 


3.79 


± 


0.12 


0.60 


± 


0.07 


-0.25 ±0.03 




ULAS1335 


18.70 


± 


0.05 


17.90 


± 


0.03 


18.15 


± 


0.12 


18.47 


± 


0.20 














0.80 


± 


0.06 


-0.25 ±0.13 


-0.32 ±0.23 




18.81 


± 


0.04 


17.90 


± 


0.01 


18.25 


± 


0.01 


18.28 


± 


0.03 


22.04 


± 


0.10 


4.14 


± 


0.10 


0.91 


± 


0.04 


-0.35 ±0.01 


-0.03 ± 0.03 



Table 2. The availa ble photometry f or each of the new ultracool dwarfs. In each case the first row gives survey photometry from UKIDSS 
LAS (apermag3; see lDve et al,| [2006y All three objects were undetected in SDSS DR6. The second row for each object gives the results 
of the follow-up photometry described in Section 12.21 



2.3 Spectroscopic confirmation 

Spectroscopy in the JHK-han&s was obtained for 
ULAS1238 and ULAS1335 us ing the Near InfraR ed Im- 
ager and Spectrometer (NIRI; [rlodapp et al.ll2003T ) on the 
Gemini North Telescope on Mauna Kea under programs 
GN-2007B-Q-26 and GN-2008A -Q-15. The InfraRed C am- 
era and Spectrograph (IRCS; iKobavashi et all |2000| ) on 
the Subaru telescope on Mauna Kea was used to obtain 
the JH spectrum for ULAS1017 and this was then fol- 
lowed with spectroscopy in the H and 7f -bands with Gem- 
ini/NIRI. Additionally, we obtained a F-band spectrum 
for UL AS 1335 usi ng the Infrared Spectr ometer And Array 
Camera (ISAAC; iMoorwood et al.l(l998T ) mounted on VLT 
UT1 and ESO, Paranal in Director's Discretionary time 
(program ID: 280.C-5067(A)). 

All observations were made up of a set of sub-exposures 
in an ABBA jitter pattern to facilitate effective background 
subtraction, with a slit width of 1 arcsec. The length of the 
A-B jitter was 10 arcsecs. The observations are summarised 
in Table [3] The NIRI observations were reduced using stan- 
dard IRAF Gemini packages. The Subaru IRCS JH spec- 
trum was also extracted using standard IRAF packages. The 
AB pairs were subtracted using generic IRAF tools, and me- 
dian stacked. In the case of IRCS, the data were found to be 
sufficiently uniform in the spatial axis for flat-fielding to be 
neglected. The ISAAC data were reduced and extracted us- 
ing the same steps, but implemented with the ESO ISAAC 
pipeline version 5.7.0. 

For all three sets of data, a comparison argon arc frame 
was used to obtain a dispersion solution, which was then ap- 



plied to the pixel coordinates in the dispersion direction on 
the images. The resulting wavelength-calibrated subtracted 
pairs had a low-level of residual sky emission removed by fit- 
ting and subtracting this emission with a set of polynomial 
functions fit to each pixel row perpendicular to the disper- 
sion direction, and considering pixel data on either side of 
the target spectrum only. The spectra were then extracted 
using a linear aperture, and cosmic rays and bad pixels re- 
moved using a sigma-clipping algorithm. 

Telluric correction was achieved by dividing each ex- 
tracted target spectrum by that of an early A or F type 
star observed just before or after the target and at a simi- 
lar airmass. Prior to division, hydrogen lines were removed 
from the standard star spectrum by interpolating the stellar 
continuum. Relative flux calibration was then achieved by 
multiplying through by a blackbody spectrum of the appro- 
priate T e ff. Data obtained for the same spectral regions on 
different nights were co-added after relative flux calibration, 
each weighted by their exposure time. 

The spectra were then joined together using the mea- 
sured near-infrared photometry to place the spectra on an 
absolute flux scale, and rebinned by a factor of three to in- 
crease the signal-to-noise, whilst avoiding under-sampling 
for the spectral resolution. The IRCS JH spectrum of 
ULAS1017 was joined to the NIRI H and J^-band spectra 
in a similar way, but the NIRI spectra were rebinned by a 
factor of three prior to joining. The resultant YJHK and 
JHK spectra for ULAS1017 and ULAS1335 are shown in 
Figure [T] where the combined spectra have been normalised 
to unity at \.27\ixa and each offset for clarity. In the case of 
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Object 


UT Date 


Integration time 


Instrument 


Spectral region 


TIT,AS1fl17 

KJ i_J JT. iJ J. \J J_ ( 


2008 T arl 2^ 


1 2x300s 


IRCS 


JH 




9008 A/Tar 1 d 




NIRI 






2008 May 11 


12x300s 


NIRI 


H 


ULAS1238 


2007 Dec 21 


4x240s 


NIRI 


J 




2008 Jan 23 


12x300s 


NIRI 


J 




2008 Jan 23 


12x300s 


NIRI 


H 




2008 Feb 18 


8x225s 


NIRI 


K 




2008 Feb 26 


6x225s 


NIRI 


K 


ULAS1335 


2007 Dec 21 


4x240s 


NIRI 


J 




2008 Jan 22 


12x300s 


NIRI 


J 




2008 Jan 22 


12x300s 


NIRI 


H 




2008 Feb 28 


20x225s 


NIRI 


K 




2008 Mar 19 


12x300s 


ISAAC 


Y 



Table 3. Summary of the near-infrared spectroscopic observations. 



I ULAS1238 




/ 1 / \ i 



1.0 1.2 1.4 1.6 1.8 2.0 2.2 2.4 

Wavelength / /j,m 

Figure 1. The NIRI JHK spectrum for ULAS1335 (bottom), 
the JH spectrum for ULAS1238 (middle). The IRCS J and NIRI 
HK spectrum of ULAS1017 is shown on the top row. 




Figure 2. The NIRI if-band spectrum for ULAS1238. The error 
spectrum is shown offset by -0.2. 



ULAS1238, we have so far been unable to obtain a _R"-band 
magnitude, so only the JH spectrum is shown in Figure [TJ 
the A'-band spectrum is shown separately in Figure [2] nor- 
malised to unity at 2.08(xm. 



2.4 Spitzer-IRAC follow-up of ULAS1335 

We obtained IRAC four-channel (3.55, 4.49, 5.73 and 
7.87 (im) photometry of the source on 2008 March 5. The 
data were obtained as part of the Spitzer Space Telescope 
Cycle 4 program #40449. All four channels have 256 x 256- 
pixel detectors with a pixel size of l'.'2xl'.'2, yielding a 
5^2x5.'2 field of view. Two adjacent fields were imaged in 
pairs (channels 1 and 3; channels 2 and 4) using dichroic 
beam splitters. The telescope was then nodded to image a 
target in all four channels(j We used exposure times of 30 s 



and a 9-position medium-sized (52 pixels) dither pattern, for 



an exposure of 270s for each channel, and a total observing 
time of 15.8 minutes to acquire all channels, slew, and settle. 

The data were reduced using the post-basic-calibration 
data mosaics generated by version 17.0.4 of the IRAC 
pipeline^ The mosaics were flat-fielded and flux-calibrated 
using super-flats and global primary and secondary stan- 
dards observed by Spitzer. We performed aperture photom- 
etry using an aperture with a 2-pixel (or 2'.'4) radius, to 
maximise the signal-to-noise ratio. The object was rela- 
tively isolated with the closest source ~ 7 arcsec away. The 
object appears single in the IRAC images, as in the case 
of the near-infrared images. To convert to total fluxes, we 
applied aperture corrections as described in Chapter 5 of 
the IRAC Data Handbook 3 of 1.205, 1.221, 1.363 and 1.571 
to channels 1 through 4, respectively. The photometry was 
converted from milli-Janskys to magnitudes on the Vega sys- 
tem using the zero-magnitude fluxes given in the IRAC Data 
Handbook (280.9, 179.7, 115.0, and 64.1 Jy for channels 1 to 
4, respectively). Photometric errors were derived from the 



1 For more information about IRAC, see 

Fazi o et al.l ((2004) and the IRAC Users Manual at 2 Information about the IRAC pipeline and data products can 

http://ssc.spitzer.caltech.edu/irac/descrip.html be found at http://ssc.spitzer.caltech.edu/irac/dh/ 
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Band 


Magnitude (Vega) 


(am 




3.55 


15.96 ±0.02 


4.49 


13.91 ±0.01 


5.73 


14.34 ±0.04 


7.87 


13.37 ±0.07 



Table 4. Mid-infrared photometry for ULAS1335. Random un- 
certainties are quoted. An additional 3% uncertainty should be 
added in quadrature to the the quoted random uncertainties to 
account for systematics. 

uncertainty images that are provided with the post-basic- 
calibration data. The magnitudes and errors are given in 
TableUJ Note that in addition to the errors in the table, there 
are absolute calibration uncertainties of 2-3 per cent. There 
are also systematic uncertainties intr oduced by pipeline de - 
pendencies of comparable magnitude (|Leggett et al . 2007b). 
We adopt the total photometric uncertainty to be the sum 
in quadrature of the values given in Table [3] plus 3%. 



do not feel that this identification is certain. Pending fur- 
ther confirmation of this feature, we will refer to this index 
as "NH3"-Jj" . Whether the a bsorp tion feature attributed to 
ammonia bv lDelorme et ah (2008) is indeed the first defin- 
ing feature of the Y class will depend on how this feature 
develops towards lower temperatures and as more objects in 
this regime are identified. 

We calculate the Wj and " NH3" — H s pectral indices 
as put forward by I Warren et al] ([2007) and iDelorme et all 
(2003) respectively. The Wj index traces the decreasing 
width of the J-band peak through the ratio of flux in a 
section of the blue slope of the peak to that in the crown. 
As such the value Wj index is driven by the strength of 
the J-band water absorption. These indices and the flux ra- 
tios from which they are derived are given in Table [5] for the 
three objects we have introduced, in addition to those for the 
T8 standard 2MASS0415 and the two previously discovered 
T8± objects. The final two columns of Table [S]give the im- 
plied type from the template comparison and the spectral 
type we adopt for each object, as discussed in more detail 
in the following sections. 



3 SPECTRAL TYPES 

Two methods for establishing spectral type are used. The 
first is by comparison of our spectra to the spectra of T spec - 
tral type standards as defined by iBurgasser et al.l |2006). 
Since the latest spectral type defined by that system is 
T8, for which the standard is 2MASS J04151954-0935066 
(hereafter 2MASS0415), we also compare our spectra to the 
two previously identified T8± objects, ULAS0034 and CF- 
BDS0059. 

The Y and K band morphologies are more strongly 
dependent on metallicity and gravity than the J and H- 
bands. This is due to the influence of pressure dependent 
collisionally induced H2 absorption in the K band (e.g. 
IBurgasser et alj|2002ft. and the red side of the Y band flux 
peak ( Allard et al.l l2008). In addition, the blue side of the 
Y band peak is shaped by the red wing of the broad KI ab- 
sorption ce ntered at 7700 A, wh ich is strongly dependent on 
metallicity (|Burrows et al.|[2002t ). As such we base our spec- 
tral types on the comparisons of the J and H band spectra 
of our objects. In this case we normalise our J and H-band 
spectra to unity at 1.27u.m and 1.58(i.m respectively, to facil- 
itate comparison of their relative shapes, rather than heights 
(which will be discussed separately). 

The second method for establishing type is the use of 
the T dwarf spectr al indices, calculated a s flux ratios, which 
are also defined by IBurgasser et all l|2006l ). Aware that some 
spectral indices may become saturated for types later than 
T8, we place our emphasis for determining type on the tem- 
plate comparison, and use this to assess the validity of the 
defi ned spectral ind i ces. 

Delormc et all (|2008h argue that ULAS0034 and CF- 
BDS0059 may represent the transition to Y spectral class 
due to the onset of broad NH3 absorption in the blue side of 
the H-band flux peak. They put forward the NH3-if ratio 
as a spectral index, and suggest that it traces the strength of 
ammonia absorption (see Table [SJ. However, given the ab- 
sence of calculated ammonia opacities in the near-infrared, 
and possible influence of nearby water absorption band we 



3.1 ULAS1017 

The spectrum of ULAS1017 has a peculiar property. It ap- 
pears to have a different spectral type in the H and K- 
band flux peaks as compared to the J-band peak. Fig- 
ure |3] shows the comparison of our IRCS±NIRI spectrum 
for ULAS1017 with the T8 standard 2MASS0415 and the T6 
standard SDSS J162414.37±002915.6 (hereafter SDSS1624). 
ULAS1017 traces the form of the T8 spectrum well in the J- 
band, but the i/-band peak matches the T6 spectrum best, 
particularly towards the the red end of this region. This is 
reflected in the calculated spectral indices for this object (see 
Table [5}. We adopt therefore the type T8p for this object. 

We do not think such morphology could arise as a re- 
sult of a T8 ± T6, or similar, unresolved binary, as flux 
from the brighter T6 would tend to fill-in the deep absorp- 
tion feature in the J-band, where we see very little flux in 
this case. Furthermore, spectral synthesis using various bi- 
nary c omponents, following the method of IBurgasser et all 
(2008), is not able to reproduce the spectral ratios seen for 
this object. The problem arises because earlier type objects, 
which can reproduce the H and K bands tend to dominate 
the light, and drive the ratios in the J band to earlier val- 
ues also. Althou gh J-band flux reve rsals have been observed 
elsewhere (e.g., lLooper et al.1120081 ). the late spectral type 
of this object would seem to preclude that interpretation. 
iPinfield et ail l|2008h identified ULAS J1150±0949, a T6.5p 
with a similarly peculiar early i/-band type (of T3 in this 
case), so such morphology may be a generic trait for some 
objects. 
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Figure 3. The IRCS+NIRI JHK spectrum of ULAS1017 (black line) overplotted with the spectra of the T8 standard 2MASS0415 and 
the T6 standard SDSS1624. The error spectrum, offset by -0.2 in the y-axis, is plotted as a black line. 
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Figure 4. The J and Jf-band flux peaks of ULAS1238 (black line), ULAS0034 (T9, sec text; green line) and 2MASS0415 (T8; red line). 
The black line that is plotted offset by -0.2 in the y-axis is the error spectrum for ULAS1238. Also indicated are the spectral regions 
integrated to form the numerator of the spectral indices discussed in the text. The region that forms the denominator is indicated with 
the letter "D". 
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Figure 5. The J and ff-band flux peaks of ULAS1335 (black line), ULAS0034 (T9, see text; green line) and CFBDS0059 (T9, see text; 
red line). The black line that is plotted offset by -0.2 in the y-axis is the error spectrum for ULAS1335. Also indicated are the spectral 
regions integrated to form the numerator of the spectral indices discussed in the text. The region that forms the denominator is indicated 
with the letter "D" . 
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J 1.14 


pi. 34 

/(A)dA 
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/(A)dA 
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pi. 60 

/(A)dA 
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pi. 60 

/(A)dA 

J 1.57 


p2.l2 

/(A)dA 

J2.O8 


J 


•1.29 

/(A)dA 

1.25 






2MASS0415 


0.030 (TS) 


0.168 (T8) 


0.31 


0.173 (T8) 


0.105 (T8) 


0.625 ± 0.003 


0.050 (>T7) 




0.134 




T8 


ULAS0034 


0.012 ± 0.006 (>T8) 


0.144 ± 0.009 (>T8) 


0.262 ± 0.006 


0.133 ± 0.010 (>T8) 


0.096 ± 0.006 (T8) 


0.516 ± 0.008 


0.091 ± 0.015 (>T7) 





.128 ± 0.003 


>T8 


T9 


CFBDS0059 


0.029 ± 0.005 (TS) 


0.164 ± 0.008 (T8) 


0.257 ± 0.004 


0.119 ± 0.008 (>T8) 


0.084 ± 0.002 (T8) 


0.526 ± 0.005 


0.128 ± 0.037 (>T7) 





.101 ± 0.002 


>T8 


TO 


ULAS1017 


0.054 ± 0.007 (T8) 


0.156 ± 0.006 (T8) 


0.347 ± 0.005 


0.274 ± 0.016 (T6) 


0.309 ± 0.011 (T6) 


0.688 ± 0.017 


0.174 ± 0.018 (T6) 





.119 ± 0.002 


T8/T6 


T8p 


ULAS1238 


0.043 ± 0.008 (T8) 


0.218 ± 0.009 (T7) 


0.284 ± 0.006 


0.145 ± 0.013 (>T8) 


0.085 ± 0.009 (T8) 


0.601 ± 0.013 


0.083 ± 0.012 (>T7) 






T8 


T8.5 


ULAS1335 


0.034 ± 0.004 (T8) 


0.179 ± 0.004 (T8) 


0.254 ± 0.003 


0.114 ± 0.004 (>T8) 


0.077 ± 0.003 (T8) 


0.564 ± 0.002 


0.094 ± 0.007 (>T7) 





.134 ± 0.001 


>T8 


TO 



Table 5. The near-infrared spectral index ratios as defined by iBurgasser et~a l. (2006) f or the three objects presented in this pap er, and for the T8 st andard 2MASS0415 and the 
previously discovered T8+ dwarfs. Also given is th e Wj index as defined by IWarrenetaL 

I fcOOTt) , and the "NH^'-H index defined bv lDelorme et aT 1 11200811 . The final two columns give 
the spectral type from template comparison in the[Burgasscr et al. (2006) scheme, and the final adopted spectral type. 
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Figure 6. The full YJHK comparison of ULAS1335 to 
ULAS0034, CFBDS0059 and 2MASS0415. 

3.2 ULAS1238 

Figure [4] compares the J and //-band flux peaks for 
ULAS1238 to those for 2MASS0415 and ULAS0034. Initial 
examination of these plots, and the spectral indices given in 
TableEJ suggest a spectral type of T8 for ULAS1238. Closer 
examination of the //-band peak reveals some similar nar- 
row features in the spectra of ULAS0034 and ULAS1238. 
Additionally, the width of the J-band peak, and the values 
of the //-band spectral indices imply that ULAS1238 is in 
fact intermediate between ULAS0034 and 2MASS0415. If 
we adopt the type T9 for ULAS0034 (see Section [373j|, then 
ULAS1238 could reasonably be classified as T8.5. However, 
the inherently small changes in near-infrared spectral mor- 
phology in this regime make such distinctions non-critical in 
the absence of additional mid-infrared spectral coverage. 

3.3 ULAS1335 

An initial inspection of the shapes of the J and //-band 
flux peaks of UL AS 1335 (see Figure [SJ indicates that it is 
of the same spectral type as ULAS0034, and CFBDS0059. 
In particular it should be noted that the strength of t he ab - 
sorption blueward of 1.565 u.m that iDelorme et alj (|2008l ) 
attribute to NH3 is similar for all three objects. Given the 
close similarity of these objects to the rest of the T dwarf 
sequence, it is appropriate that these objects form an ex- 
tension of that sequence, and we assign to them the sub- 
type T9. Figure [6] shows the full YJHK comparison of 
ULAS1335, CFBDS0059, ULAS0034 and 2MASS0415. Al- 
though the shapes of the //-band peaks are similar for the 
three objects, their height relative to the J-band peak varies 
significantly. Specifically, ULAS1335 and CFBDS0059 both 
exhibit much stronger //-band peaks than either ULAS0034 
or 2MASS0415. 

In Figure [7] we plot the J and H band spectral indices 
against the spectral types derived by template comparison 
for objects with types T6 to T9 (see below). It is apparent 
that the water and methane J-band indices previously used 
for the classifying T dwarfs are degenerate with type for the 
subtypes T8 and T9. Examination of Figure [7] indicates that 



the same may be true for the CH4-ZZ index, although the 
H2O-Z/ index is still able to distinguish these objects. The 
new indices "NH3"-// and Wj both appear able to distin- 
guish the T9 dwarfs from the earlier type objects, however 
the "NH3"-// index is degenerate with type for the earlier 
objects. 

In Figure [8] we plot the "NH3"-// inde x again st the Wj 
spectral index suggested by IWarren et all (|2007l ). It is ap- 
parent that the latest objects are well separated from ^T8 
objects in this plot. However, the values of Wj and "NH3"- 
H do not show a consistent correlation for the latest objects. 
Figure U] shows the H2O-// index plotted with Wj. Again, 
the three latest objects are well separated from the other 
objects, however in this case the values of the two indices 
show a clear trend (with the exception of the T8p dwarf), 
both decreasing together. The scatter in "NH3"-// values 
in Figure [S] compared to the correlation shown in Figure 
thus suggests that although "NH3"-// is useful for identify- 
ing T8+ objects, it is not an effective typing index for very 
late objects in the range we are considering, and does not 
conveniently strap on to the previous T dwarf sequence. 

We suggest that the appropriate spectral indices for the 
extension of the T-spectral class to T9 are H2O-Z/ and W.j. 
To provide continuity with the earlier T-indices we define the 
range of Wj and H2O-Z/ fo r types T7-T9, matchin g our val- 
ues to the H2O-Z/ values of lBurgasser et ahl (|2006h for types 
T7 and T8, and using the apparent correlation in Figure[9]to 
arrive at appropriate Wj values. Our proposed set of spec- 
tral indices for the latter portion of the T spectral sequence 
are summarised in Table [6] In this scheme, ULAS1335, 
as well as ULAS0034 and CFBDS0059, are all T9 dwarfs, 
which are our adopted types. Since ULAS1335 is signifi- 
cantly brighter than both ULAS0034 and CFBDS0059, we 
also propose ULAS1335 as the template for the T9 subtype. 

In Section r5.2l we esti mate a r e g for UL AS 1335 which is 
marginally cooler than the I Warren et all (120071 ) estimate for 
ULAS0034. This temperature sequence is well represented 
by the locations of the objects in Figure [5] but not by the 
sequence shown in Figure(8] If the "NH3"-// index is indeed 
tracing a broad ammonia absorption band, this suggests two 
things. Firstly, it appears to develop at odds with the water 
absorption in the J-band (which the Wj feature traces) and 
the H2O-// absorption feature, Secondly, if the relative tem- 
peratures of ULAS0034 and UL AS 1335 are to be believed, 
it does not correlate well with temperature. 

There is some theoretical justification for expecting 
near-infrared ammonia absorption to be only weakly depen- 
dent on T e ff. Ammonia chemistry is governed by upmixing, 
and in practice the mixing ratio depends most strongly on 
the ammonia abundance at the quench level. Since this is 
at a nearly fixed temperature, ~2000 K, the remaining vari- 
ables that will affect it are pressure and composition. Thus, 
we expect the ammonia mixing ratio, and by extension the 
strength of ammonia absorption in the near-infrared, to be 
far more sensitive to gravity and metallicity than T a g . Given 
this behaviour, there is exciting potential for the "NH3"-// 
index as an additional diagnostic tool for very cool T dwarfs, 
and its development will be followed with interest as more 
more objects in the very low-temperature regime are iden- 
tified. 



10 Ben Burningham et al. 




T subtype T subtype T subtype 

Figure 7. Spectral index versus T-subtype for T6-T9 dwarfs. Asterisks indicate "normal" dwarfs, whilst the diamonds indicate two T8p 
dwarf s: ULAS1017 and 2MA SS J 07290002-3954043. Index values for these objects are drawn from iBurgasser et al.l l|2006l)IWarren et al] 
(|2007l l. lDelorme et al.lll2008ll and lLooper et alJfeOOTft . or are calculated from the objects' spectra supplied by these authors. Uncertainties 
are smaller than the symbol sizes. 



NIR spectral type 


H2O-J 


CH4-J 


Wj 


H2O-H 


CH4--H" 


CH4-K 


Standard 


T7 


0.07-0.13 


0.21-0.28 


0.35-0.40 


0.20-0.26 


0.15-0.25 


<0.13 


2MASS0727 


T8 


<0.07 


<0.21 


0.28-0.35 


0.14-0.20 


<0.15 




2MASS0415 


T9 






< 0.28 


< 0.14 






ULAS1335 



Table 6. Our proposed spectral indices for the latter-most portion of the T spectral sequence. Also given are the spectral type standards 
for each type. 



4 SEARCHING FOR COMMON PROPER 
MOTION COMPANIONS 

It has been known for some time that the degree of multi- 
plicity amongst very young stars is gr eater than that of the 
more evolved field star populations (|Duauennov fc Mavorl 
ll99ll ; lLeinert et al.lll993h . and thus that the majority of bi- 
nary systems form together in their nascent clouds. Binary 
components can therefore generally be assumed to share the 
same age and composition, and wide common proper motion 
companions to our new objects could thus provide valuable 
constraints if they are to be found. In order to address this 
possibility and more generally assess the kinematics of our 
objects, we calculated proper motions for our targets using 
the highest signal-to-noise multi-epoch observations avail- 
able. In each case this was the WFCAM and UFTI J-band 
images. 

The IRAF routines geomap and geoxytran were used to 
transform between the available multi-epoch images, using 
an average of 12 reference stars. This allowed any motion of 



the T dwarfs to be accurately measured. Proper motion un- 
certainties were estimated from centroiding accuracies com- 
bined with the residuals associated with our derived trans- 
formations. The calculated proper motions are given in Ta- 
ble [7] We also estimate Vtan based on our distance estimates 
derived in Sectio n [7] These e s timat es lie within the range of 
Vtan observed by IVrba et al.l (|2004T l for L and T dwarfs. 

We independently searched the Supercosmos Sky Sur- 
vey database, and the Simbad database (accessing the Hip- 
parcos and Tycho catalogues amongst others), for common 
proper motion companions (within ler) of the three new T 
dwarfs, thus searching for wide binary companions over a 
wide range of mass. The search was made within a projected 
distance of 20,000 AU from each target, calculated assuming 
the minimal distance within our derived constraints (see Sec- 
tion [TJl to ensure a spatially complete search. In the case of 
ULAS1238 we used a lower distance estimate of 20 pc based 
on comparison with ULAS0034. SDSS and UKIDSS pho- 
tometry were then retrieved for any common proper motion 
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Figure 8. The "NH3"— H index versus Wj index. Uncertainties 
are not displayed as they are smaller than the symbol size. 
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Figure 9. The H2O-H index versus Wj index. Uncertainties are 
not displayed as they are smaller than the symbol size. 



objects, and from these a spectral type range was estimated, 
and an approximate distance modulus inferred. Sources with 
distance moduli that were inconsistent with the inferred dis- 
tances to our newly discovered dwarfs (see Section [7} could 
then be rejected as potential companions. Our search pro- 
duced no potential companions for any of our targets. 



5 ATMOSPHERIC PARAMETERS 

In the following two sections we estimate the atmospheric 
parameters of ULAS1017 and ULAS1335 through model 
comparisons, and, in the case of ULAS1335, consideration 
of the near-infrared to mid-infrared spectral energy distri- 
bution (SED). The model spectra we employ are a recent 
realisation of the BT-Settl models, for whic h detailed de- 
script ions of the input p hysic s are given in IWarren et alj 
i|2007h and lAllard et al.l (|2008h . There are uncertainties in 
these models that limit their use significantly for simplistic 
comparisons, but by limiting our analysis to relative com- 
parisons over key spectral ranges, and some extrapolation 
of theoretical trends, we are able to derive a set of "best- 
guess" properties based on the existing theoretical and ob- 
servational data. 

We do not perform a detailed analysis for UL AS 1238 
due to the absence of Jf-band photometry that precludes 
flux calibration for this part of the spectrum. Instead we 
speculate that its properties are likely intermediate between 
2MASS0415 and ULAS0034 since its spectral type is also 
intermediate. However, we do not rule out the possibility of 
a considerably lower T c s (see Section O. 



5.1 ULAS1017 

To estimate the atmospheric parameters for ULAS1017 we 
repeat the (Wj,K/J) analysis described by IWarren et alj 
i200H\ for ULAS00 34, and applied to CFBDS0059 by 
iDelorme et~ai1 (120081 ). To this end we have calculated Wj 
and Kj J indices for a grid of solar metallicity BT-Settl mod- 
els with 4.20 < log g < 5.50, and 600.fr < T aS < 900K. 
We have anchored the resulting (Wj,K/ J) grid such that 
T c ff = 750K and log g — 5.0 for the empirical (Wj,K/J) co- 
ordina tes of 2MASS0415, c onsistent with the detailed anal- 
ysis of lSaumon et all l|2007l ). Placing ULAS1017 on this grid 
implies a T e g = 800K and log g=5.5. 

In Figure [TT] we compare the spectrum of ULAS1017 
to model spectra that bracket its location in the (Wj,K/ J) 
plane. It appears that the model spectra generally fit the 
spectrum of ULAS1017 well. We note, however, that it is 
the model for T e g = 800K and log g = 5.0 that provides the 
best fitting spectrum, in contrast to the (Wj,K/J) analysis 
that suggests log g = 5.5. 

The location of ULAS1017 in Figure QJj] is very 
close to that of Gl 570D. It is worth noting that us- 
ing de tailed analysis with c onsideration of non-equilibrium 
effects lSaumon et al.l (|2006l ) derive parameters for Gl 570D 
of T efl = 800 - 820K and log g = 5.09-5.23. This is very 
similar to our best value from spectral comparison for 
ULAS1017, but similarly at odds with the parameters from 
simple (Wj,K/J) analysis. 

Note also that this form of analysis can only provide ap- 
proximate estimates of the atmospheric parameters, and it 
must be remembered that varying estimates of surface grav- 
ity will be degenerate with respect to varying metallicity. 
The lack of Y-band coverage in our spectrum of ULAS1017 
prevents us from making any estimate of its metallicity rel- 
ative to other objects based on the form of its F-band peak. 
Furthermore, the form of the i/-band flux peak, resembling 
that of a T6 dwarf adds an additional element of uncertainty 
to our "best-guess" properties for this object. 
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Object 


1st Epoch Coordinates 


Epoch 1 


Epoch 2 


^acosS 




Vtan 




a 


S 






mas yr 


mas yr — 1 


kms — 1 


ULAS1017 


10 17 21.40 


+01 18 17.9 


2007-03-17 


2008-01-17 


110 + 80 


40 + 80 


4-58 


ULAS1238 


12 38 28.51 


+09 53 51.3 


2007-02-13 


2008-01-25 


418 ± 40 


21 + 40 




ULAS1335 


13 35 53.45 


+11 30 05.2 


2007-04-21 


2008-01-16 


160 ± 50 


-290 ± 50 


10-21 



Table 7. Proper motion estimates for the three new T dwarfs. Uncertainties reflect the residuals in the transformations between the two 
images in each case. Also given are estimates for Vtan for ULAS1335 and ULAS1017 (see Section [7] for distances). 
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Figure 10. Wj versus K/J indices for a grid of solar metallicity 
BT-Settl model spectra. The grid has been normalised by match- 
ing the T c fj = 750.R", log g = 5.0 model indices to the empirical 
index values for 2MASS0415, which lies at (1.0, 1.0). 



5.2 ULAS1335 

We will begin by placing ULAS1335 on the (Wj,K/J) grid 
(see Figure [TUJ , which implies T e ff « 650K and log g ~ 4.5. 
To assess the credence of these theoretical predictions, Fig- 
ure \T%\ shows the comparison of UL AS 1335 with the BT- 
Settl synthetic spectra for the atmospheric parameters that 
bracket its location in the (Wj,K/J) plane for solar metal- 
licity. It is clear that these synthetic spectra are not a good 
fit to the data. The most obvious discrepancy is the height 
of the //-band flux peak. It is widely acknowledged that the 
shape of the //-band peak is poorly reproduced by model 
spectra, a fact that is attributed to incomplete methane line 
lists. However, in the cases of ULAS0034, CFBDS0059 and 
2MASS0415, the observed height of the //-band flux peak is 
well represented by the model spectra for the atmospheric 
paramet ers estimated from their location in the ( Wj,K/J) 
plot (see I Warren et al. l l2007l ; lDelorme" et al. I [20081 . figures 6 
and 7 respectively) It is thus clear that despite many sim- 
ilarities between ULAS1335 and the other two T9 dwarfs, 
its //-band peak is, by comparison, significantly brighter. 
Some of the properties of UL AS 1335 are thus presumably 
somewhat different to the other two T9 dwarfs. 
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Figure 11. A comparison of model spectra for the parame- 
ters bracketing the location of ULAS1017 in the solar metallicity 
(Wj,K/J) plane. 



We can make a first assessment of possible differences 
in metallicity by considering the shape of the Y-band peak. 
All available theoretical models indicate that the shape of 
this flux peak will be sensitive to [M/H], although there is 
some ambiguity as to exactly how compositional changes 
will affect its shape. A comparison of the Y-band peaks of 
the three T9 dwarfs in Figure [6] reveals them to have an 
extremely similar Y-band spectral morphology, and it thus 
seems unlikely that these objects have significantly different 
metallicity. 

Figure [T2l shows what changes we might expect for vary- 
ing log g and T c g . In the bottom plot it can be seen that 
although a higher gravity can increase the brightness of the 
//-band peak relative to the J— and Y-band peaks, it has 
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the opposite effect on the Tf-band flux peak. These trends 
thus indicate that it is not possible to theoretically account 
for the relatively bright .ff-band peak of ULAS1335 com- 
pared to the other T9 dwarfs, simply by varying log g alone. 

We can obtain a better match between theoretical and 
observed spectra if we adopt a lower T e fr for ULAS1335. The 
top plot in Figure [12] demonstrates the trend in the BT- 
Settl model spectra, and shows that the relative brightness 
of the if-band peak increases compared to the Y- and J- 
band peaks as T c g decreases below 700K. However, note 
that although there is some change to the relative strength 
of the A'-band peak, it is not as drastic as when one varies 
log g. The overall near infrared morphology of ULAS1335 is 
thus best explained (through current theory) if this object 
has a somewhat lower T e g than ULAS0034. 

In Figure [13] we compare the spectra of UL AS 1335 
to a range of model spectra with different metallicities 
and gravities for our lowest available temperature models, 
Toft = 600K. In the top plot, the combination log g = 5.0, 
T c ft = 600K model spectra all estimate the height of the H- 
band peak well. The sub-solar metallicity model is a poor 
match everywhere else however. The high-metallicity model 
gives a reasonable approximation to the observed spectrum, 
but as expected it under-estimates the K-banA flux some- 
what, and has a significantly wider and blue- ward shifted 
Y-band peak. 

In the bottom plot, the log g =4.25-4.50 model spectra 
provide a reasonable fit to the whole near infrared spectra. 
However, whilst reproducing the Y-, J-, and A'-band peaks, 
these spectra still somewhat under-estimate the brightness 
of the if-band peak. Although the low T e ff of this model 
has gone some way to addressing the relative brightness of 
this peak, the model trends suggest that we should obtain a 
better fit by going to even lower T e ff. 

Investigation of the near-infrared to mid-infrared SED 
for ULAS1335 reveals that is extremely red. Figure lFfl shows 
the l-8u.m SEDs for ULAS1335, ULAS0034, 2MASS0415 
and Gl 570D. Whilst the form of the SEDs are simi- 
lar, ULAS1335 is significantly redder than ULAS0034 over 
the 2 - 4.49 u.m range. ULAS1335 has H - [4.49] = 
4.34 ± 0.04, the red dest value yet observed for a T dwarf. 
IWarren et ail |2007l ) identified a correlation between the H— 
[4.49] colour and Teff for T dwarfs with well determined 
T e ff i Golimowski et ail [2004). In fact the correlation for T 
dwarfs with T e ff = 800 — 1250K is particularly tight, with 
an RMS about the best-fitting linear relation of ±28K. In 
Figure [TSI we show this linear fit, and indicate the _H-[4.49] 
colours of 2MASS0415, ULAS0034 and ULAS1335. Theo- 
retical _tra£ks_arealso_over-plotted for a recent realisation 
of the lMarlev et all (l2002f) model at mospheres for the range 
500-800K (|Marlev fc Saumonl 12008) 1 . Tracks are shown for 
4.48 < log 5 < 5.48 and -0.3 < [M/H] < +0.3, bracketing 
the expected properties of ULAS1335. We have anchored the 
models for log g = 5.0 and solar metallicity to the colours 
of the coolest T dwarfs with well determined temperatures. 
It can be seen that the theoretical models are reasonably 
supportive of an extrapolation of the relation between H- 
[4.49] and T e g for such log<? and [M/H] ranges, and this 
extrapolation suggests that ULAS1335 has T e « ~550-600K. 

In the case of ULAS0034, the temperature inferred by 
(Wj,K/ J) analysis was reasonably consistent with that in- 
ferred from an extrapolation of the T e ft—H— [4.49] relation. 



UWS1335 

BT Settl 700K [M/H] = 0.0 log(g) = 4.5 
BT Settl 650K [M/H] = log(q) = 4.5 
BT Settl 600K [M/H] = 0.0 log(g) = 4.5 
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BT Settl 650K [M/H] = 0.0 log(g) = 4.25 
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Figure 12. A comparison of model spectra for parameters brack- 
eting the location of ULAS1335 in the solar metallicity (Wj,K/ J) 
plane. 



However, in the case of ULAS1335 we have seen that the 
two estimates differ by ~ 100K. The poor agreement be- 
tween the model spectra for 650 K and the observed spec- 
trum of ULAS1335 cast doubt on the temperature derived 
from the (Wj,K/ J) method. It is likely that this disagree- 
ment is principally due to the fact that T9 dwarfs are in the 
temperature range near the limit of where the near-infrared 
model spectra are valid. Since no model spectra cooler than 
600 K are plotted, a degeneracy in (Wj ,K/ J) for these tem- 
peratures cannot be ruled out. Indeed, Figure ITUl indicates 
that there is a predicted distortion in the grid (which will 
result in some degeneracy) for the lowest gravity and tem- 
perature combinations currently considered. 

The strong, relatively tight T g correlation offered by 
observed T dwarf 7/-[4.49] colours combined with the the- 
oretical expectations for a reasonably steady shift of flux 
from the near-infrared into the mid-infrared with decreas- 
ing Toff, lends weight to this colour as a better temperature 
indicator for such very late objects, and we adopt the tem- 
perature implied by this colour rather than that implied by 
the (Wj,K/J) method. We also suggest that the (Wj,K/J) 
method be used with caution in the T e ff regime near the 
limit of current atmospheric models. 

To summarise, our comparison to model near infrared 
spectra suggested that ULAS1335 has T c g < 600K. This 
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ULAS1335 

BT Settl 600K [M/H] = + 0.3 log(g) = 5.0 

BT Settl 600K [M/H] = log(g) = 5.00 

BT Settl 600K [M/H] = -0.5 log(g) = 5.00 
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ULAS1335 

BT Settl 600K [M/H] = 0.0 log(g) = 5.00 

BT Settl 600K [M/H] = 0.0 log(g) = 4.50 

BT Settl 600K [M/H] = 0.0 log(g) = 4.25 



Wavelength fim 

Figure 13. A comparison of model spectra with ULAS1335 for 
differing gravity and metallicity combinations at T e ff = 600K. 



is clearly supported by the the near-to-mid infrared SED, 
which suggests T e g ~ 550 — 600K from a theoretically in- 
formed extrapolation. Our "best-guess" for the temperature 
of ULAS1335 is thus T cff -550-600K. 

This temperature is stated along with an important 
caveat. We note that preliminary i n vestig ations of re- 
cent r^alis^ti£ns_of_the Marlev et al.l (|2002h model spec- 
tra (|Marlev fc Saumonll2008l ) do not show the same trends 
with varying parameters in the near-infrared as the BT-Settl 
models. Although we defer a detailed comparison of these 
model spectra with the BT-Settl models to a future work, 
we highlight those differences of particular relevance here. 
Most notable is the absence of a trend of increasing height 
of the H band peak with decreasing T e g. In addition the 
height of the H band peak appears more strongly depen- 
dent on metallicity than in the BT-Settl models, increasing 
with decreasing metallicity. 

Finally, considering the absence of mid-infrared pho- 
tometry for CFBDS0059, and the reliance on the (Wj,K/J) 
method for its temperature estimate, we also note that 
no strong conclusion regarding its temperature relative to 
ULAS1335 and ULAS0034 should be drawn. We cannot rule 
out the possibility that CFBDS0059 is cooler still. 
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Figure 14. The 1-8 urn SEDs for ULAS1335,ULAS0034 
llWarren et alj|2007t) . 2MASS0415 and Gl 570D. The photome - 
try for the latter two o bjects, taken from iLeggett et al.l (|2002h 
and iPatten et al ] l|2006l l. has been offset by +2.5 magnitudes to 
allow comparison with the T8+ dwarfs. Uncertainties are of com- 
parable size to the symbols. 
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Fi gure 15. The T e g - _ff-[4.49] relation for T dwarfs as defined 
bv lWarren et all <2007h . The shaded regions indicate the colours 
(Her) for the coolest T dwarfs with IRAC magnitudes. The 
straight red line indicates the empirical fit the colours of warmer 
T dwarfs. Also plott ed, as black lines are the o retica l tracks from 
iMarlev et all j2002f ) and iMarlev fc Saumonl l|2008l l. Values for 
logg and metallicity are indicated in brackets. 
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6 COMPARISON WITH THE JOVIAN 
SPECTRUM 

We have also examined our T dwarf spectra in more detail, 
to try and identify new fe atures that might be u seful in this 
new T e ff regime. We follow lLeggett et al.l l|2007al ) and use the 
spectrum of Jupiter along with transmission spectra of am- 
monia and methane as a basis for comparison, which we plot 
in Figure ITol alongside ULAS0034, CFBDS0059,ULAS1238 
and ULAS1335. Although the metallicity of the Jovian at- 
mosphere is much higher, and the gravity much lower, than 
that expected for our T dwarfs, this comparison with a much 
cooler atmosphere is still interesting. 

We draw attention to the narrow feature near 1.23(i.m 
(left-most vertical green line in Figure I16|l which appears 
to be shared by Jupiter, ULAS1238 and ULAS1335. This 
should not be confused with the K I doublet at ~ 1.24 — 
1.25(i.m. We attribute this feature to methane, based on 
comparison with the l abora tory transmission spectrum of 
ICruikshank fc Binderl (jl969h . We also highlight the narrow 
feature at ~ 1.272[i.m (right-most vertical green line), which 
also appears to be common between ULAS1238, ULAS1335 
and Jupiter. This may be due to ammo nia, based on com - 
parison to the transmission spectrum of llrwin et all (QUI). 
We note, however that the sharp feature blueward of this 
(indicated by a vertical red line); also likely due to ammo- 
nia is apparently absent from the spectra of the T dwarfs. 



7 PHYSICAL PARAMETERS 

Since metallicity is thought to only weakly impact the 
deriv ation of fundamenta l parameters from T c g and log g 
(e.g.. ISaumon et al.ll2007f). we use the solar metallicity rela- 
tions summarised by Burrows et al.l (|200lT > to estimate the 
mass, age and radius for ULAS1017 and ULAS1335 from 
our T e ff and log g estimates. 

For ULAS1017, the extreme cases from the model com- 
parison, T c ff = 750 — 850K and log g — 5.0-5.5 correspond to 
an age of 1.6-15 Gyr and a mass of 33-70 Mj. By estimating 
its radius from the same set of relations, and normalising the 
corresponding atmospheric models to the J-band flux peak, 
we can estimate the distance to ULAS1017 to be in the range 
31-54pc. 

For ULAS1335, we use our best estimate of 
Teff = 550 — 600K and, since we suspect low-gravity for this 
object, we adopt a conservative range of log g = 4.5 — 5.0 
for its surface gravity. These atmospheric parameters only 
weakly constrain the age to a range of 0.6-5.3 Gyr and the 
mass to 15-31 Mj. 

In the absence of model spectra for T e g < 600K, we 
have estimated the bolometric flux for ULAS1335 using a 
pseudo-synthetic spectrum. This was constructed by brack- 
eting our flux-calibrated YJHK spectrum with the model 
spectrum for T e g = 600K and log g = 4.5, scaling the syn- 
thetic parts with our z band photometry in the optical and 
our IRAC [4.49] photometry in the mid-infrared. We then 
integrated this this pseudo-synthetic spectrum over all wave- 
lengths to estimate the bolometric flux from ULAS1335. The 
uncertainty in this estimate is dominated by uncertainty in 
the photometry used to calibrate the spectrum. These un- 
certainties range from 10% for the optical photometry to 



< 5% for our near-infrared photometry. These regions con- 
tribute ~ 0.5% and ~ 25% to the bolometric flux estimate 
respectively. The mid-infrared, which contributes ~ 75% of 
the bolometric flux, has an uncertainty of 4% in its [4.49] 
magnitude, which we use to scale the synthesised spectrum. 
Allowing for additional uncertainty in the model spectrum 
match to the true mid-infrared spectrum, we estimate a to- 
tal uncertainty in our bolometric flux to be ~ 15%, and the 
bolometric flux to be F bol = 3.35 ± 0.5 x 10~ 16 Wm" 2 . 

The uncertainty in our distance estimate for UL AS 1335 
is dominated by the uncertainty in T e ft (and thus luminos- 
ity) . For Teff = 550 - 600K and log g = 4.5-5.0 we estimate 



a luminosity of —6.18 ^ log (L */L&) 



(by in- 



ferring radii using the relations in iBurrows et al.ll200lh and 
thus a distance of 8-12 pc. Table [8] summarises the inferred 
properties for the T dwarfs identified in this paper. 



8 SUMMARY 

We have identified three very late-type T dwarfs in the 
UKIDSS LAS DR3: ULAS1017, ULAS1238 and ULAS1335, 
for which we have adopted the spectral types T8p, T8.5 and 
T9 respectively. ULAS1017 has a peculiar spectrum, with a 
J-band typical of a T8 dwarf, but H and i^-band peaks of 
a T6. In assigning spectral types to the our targets, we have 
defined the extension of the T-dwarf sequence to the type 
T9, using the H 2 0-if and Wj indices, and ULAS1335 as 
the spectral standard. 

To estimate atmospheric parameters for ULAS1017 and 
ULAS1335, we have performed a detailed comparison with 
BT-Settl model spectra. For ULAS1017, we have estimated 
that its temperature lies in the range 750 ^ T e ft ^ 850, with 
log g = 5.0 — 5.5 (assuming solar metallicity). 

Our model comparison for ULAS1335 has suggested a 
Teft that is cooler than the coolest available BT-Settl mod- 
els (T c ff < 600K), and it may have low surface gravity pos- 
sibly combined with high metallicity. As a result we have 
based our "best-guess" of the atmospheric parameters on 
an examination of the near- to mid-infrared SED. We have 
found ULAS1335 to be the reddest T-dwarf yet identified, 
with H - [4.49] = 4.34 ± 0.04. By extrapolating the em- 
pirical Teff - i/-[4.49] relation to lower T ff, and guided 
by model trends, we have estimated that ULAS1335 has 
Teft ~ 550 - 600K. This estimate suggests that ULAS1335 
is the coolest brown dwarf yet discovered, although we await 
the determination of its distance by parallax with interest. 
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Figure 16. Spectra for ULAS0034, CFBDS0059, ULAS1238 and ULAS1335 compared with the near-infrared spe ctrum of Jupiter(J. T 
Rayn er, M. C. Cush ing, fc W. D. Va cca 2008, in preparation) and laboratory transmission spectra of methane (Cruikshan k fe Binder! 
Il969t) and ammonia lllrwin et al Jl99Sh . The vertical green dashed lines highlight narrow features that appear to be common to ULAS1238, 
ULAS1335 and Jupiter. The vertical red dashed line highlights an interesting feature in the Jovian spectrum that is absent from the 
T dwarf spe ctra (see text). The v ertical blue dotted line indicates the location of the red edge of the broad absorption attributed to 
ammonia bv lDelorme et al (2008). 



Object 


SpType 




log g 


mass 


age 


distance 






K 


log(cm s_ 2 ) 


Mj 


Gyr 


pc 


ULAS1017 


T8p 


750-850 


5.0-5.5 


33-70 


1.6-15 


31-54 


ULAS1335 


T9 


550-600 


<5.0 


15-31 


0.6-5.3 


8-12 



Table 8. A summary of our estimates for the properties of ULAS1017 and ULAS1335. Estimates for ULAS1238 have been generally 
neglected in the absence of full flux calibrated spectral coverage. Distances assume that the dwarfs are single objects. 
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